Low temperature broken symmetry phases of spiral antiferromagnets 
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We study Heisenberg antiferromagnets with nearest- (Ji) and third- (J3) neighbor exchange on 
the square lattice. In the limit of spin S — » 00, there is a zero temperature (T) Lifshitz point at 
J3 = \Ji, with long-range spiral spin order at T = for J3 > \Ji. We present classical Monte 
Carlo simulations and a theory for T > crossovers near the Lifshitz point: spin rotation symmetry 
is restored at any T > 0, but there is a broken lattice reflection symmetry for < T < Tc ~ 
(Ja — jJ{)S'^ . The transition at T = Tc is consistent with Ising universality. We also discuss the 
quantum phase diagram for finite S. 



Frustrated antiferromagnets have recently attracted 
much interest in connection with the possibihty of stabi- 
hzing unconventional low temperature (T) phases, pos- 
sibly falling outside the known paradigms of condensed 
matter physics 0. In this respect, the most interest- 
ing systems are those where the combined effect of frus- 
tration and quantum fluctuations is strong enough to 
prevent the onset of magnetic ordering, thus stabilizing 
'exotic' quantum disordered ground states A very 
promising candidate for such a spin-liquid phase has been 
found in the two-dimensional Ji— J3 model defined by the 
following Hamiltonian 
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where. Si are spin-S' operators on a square lattice and 
Ji, J3 > are the nearest- and third- neighbor anti- 
ferromagnetic couplings along the two coordinate axes. 
For this model, early large N computations, 01 and re- 
cent large scale Density Matrix Rcnormalization Group 
(DMRG) calculations for S* = 1/2 have suggested 
the existence of a gapped spin liquid state with ex- 
ponentially decaying spin correlations and no broken 
translation symmetry in the regime of strong frustration 
(J3/J1 - 0.5). 

This paper will describe properties of the above model 
for large S, and discuss consequences for general S. Our 
results, obtained by classical Monte Carlo simulations 
and a theory described below, are summarized in Fig^ 
for the limit S 00. There is a T = state with long- 
range spiral spin order for J3 > iJi. We establish that 
at < T < Tc ~ (J3 - I Ji)5'^ above this state there is 
a phase with broken discrete symmetry of lattice reflec- 
tions about the x and y axes, while spin rotation invari- 
ance is preserved. This phase has 'Ising nematic' order. 
We present strong numerical evidence that the transi- 
tion at Tc is indeed in the Ising universality class. Such 
Ising nematic order was originally proposed in Ref. Q 
for S* = 1/2 in a T = spin liquid phase described by a 
Z2 gauge theory . Thus the same Ising nematic order 



FIG. 1: Phase diagram of H in the limit S —> 00. The 
shaded region has a broken symmetry of lattice reflections 
about the x and y axes, leading to Ising nematic order. The 
Ising transition is at the temperature Tc ~ {J3 — jJi )S^ . The 
spin correlation length, ^spin, is finite for all T > 0, with the 
T dependencies as shown, with c/2 = c' = STrjJs — |Ji|; the 
crossovers between the different behaviors of ^spin are at the 
dashed lines at T ~ | J3 — jJi\S'^. Spin rotation symmetry is 
broken only at T = where ^spin = 00. There is no Lifshitz 
point at finite S because it is pre-empted by quantum 
effects within the dotted semi-circle: here there is a T = 
spin gap A ~ Sexp(— cS) and spin rotation symmetry is 
preserved. This semicircular region extends over T ~ IJ3 — 
jJi\S ~ A. Further details on the physics within this region 
appear at the end of the paper. 



can appear when spiral spin order is destroyed either by 
thermal fluctuations (as in the present paper: see Fig^ 
or by quantum fluctuations (as in Ref. Our large S 
results are therefore consonant with the possibilityof a 
spin liquid phase at S = 1/2 as described in Ref. Plljl; 
we will discuss the quantum finite S phase diagram fur- 
ther towards the end of the paper. We also suggest that 
discrete lattice symmetries may play a role near other 
quantum critical points with spiral order 0. 

Broken discrete symmetries have also been discussed 
0,0] in the context of the J1 — J2 model, with first- and 
second-neighbor couplings on the square lattice. How- 
ever this model has only coUinear, commensurate spin 
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FIG. 2: The two different minimum energy configurations 
with magnetic wave vectors Q = {Q,Q) and Q* = {Q,—Q) 
with Q = 27r/3, corresponding to J3/J1 = 0.5. 



correlations, and this makes both the classical and quan- 
tum theory quite different from that considered here. As 
will become clear below, the spiral order and associated 
Lifshitz point, play a central role in the structure of our 
theory and in the T dependence of observables. 

We begin by recalling ^3 the ground states of Ti at 
S ^ 00. There is conventional Neel order with magnetic 
wavevector Q = (7r,7r) for J3/J1 < \. For J3/J1 > \ 
the ground state has planar incommensurate antiferro- 
magnetic order, with a pitch vector depending on the 
frustration ratio. More precisely, the magnetic ordering 
can be described by a wavevector Q = (Q, Q) with Q de- 
creasing from TT as J3/J1 > i and approaching Q = 7r/2 
monotonically for J3/J1 — > cx); at J3 = 00 we obtain four 
decoupled Neel lattices. The spiral order is in general 
incommensurate for i < J3/J1 < 00, with the exclusion 
of J3/ Ji = 0.5 where Q = 27r/3, corresponding to an an- 
gle of 120° between spins (See Fig. |2J). Interestingly, for 
each spiral state with Q ~ {Q, Q) there is an energeti- 
cally equivalent configuration with Q* = {—Q, Q), which 
is distinct from the first one for Q ^ n. This state can- 
not be obtained from the one with wavevector Q by a 
global rotation of all the spins. Instead, the two config- 
urations are connected to each other by a global rota- 
tion combined with a reflection about the x ot y axes, so 
that the global symmetry of the classical ground state is 
0(3) XZ2, with an additional two-fold degeneracy beyond 
that of the Neel case. 

One of the main claims in Fig ^ is that the bro- 
ken Z2 symmetry survives for a finite range of T > 0, 
while continuous 0(3) symmetry is immediately restored 
at any non-zero T, as required by the Mermin- Wagner 
theorem We established this by extensive Monte 

Carlo simulation using a combination of Metropolis and 
over-relaxed algorithm for periodic clusters of size up to 
M = 120 X 120, and for several values of J3/J1 between 
0.25 and 4. Indeed, the presence of a finite T phase tran- 
sition is clearly indicated by a sha rp p eak of the specific 
heat which is illustrated in Fig.|31 |l3| This sharp feature 
is to be contrasted to the broad maximum displayed by 
the same quantity for J3/J1 < j, i.e., when the classical 
ground-state displays ordinary Neel order. In particular. 
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FIG. 3: T dependence of the specific heat for J3/J1 — 0.5. 
Different symbols refer to different clusters with liner size be- 
tween L — 24 and L = 120. Data for J3/J1 =0.1 are shown 
for comparison (full dots and dashed line). Inset: size-scaling 
of the maximum of the specific heat. 
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FIG. 4: Bottom: T dependence of the order parameter, a, 
(see Eq. ||5J) for different cluster sizes and J3/J1 = 0.5. The 
inset shows the data collapse according to the scaling hypoth- 
esis with Ising exponents /3 = 1/8 and 1^ = 1, and Tc = 0.303. 
Top: Temperature dependence of the susceptibility of a for 
J-i/Ji =0.5. The inset shows the size scaling of the maximum 
of the susceptibility. 



the maximum of the specific heat is consistent with a 
logarithmic dependence on system size (see the inset of 
Fig. 121) corresponding to a critical exponent a = 0, in 
agreement with Ising universality. 

This critical behavior can be directly related to the 
broken lattice reflection symmetry by studying an ap- 
propriate Ising ncmatic order parameter. From the sym- 
metries of Fig 121 we deduce that the order parameter is 



a = 1/M (^^ a„) with 



Si • S3 — S2 • S- 



(2) 



where a labels each plaquette of the square lattice and 
(1, 2, 3, 4) are its corners. The variables Oa are identically 
zero for a Neel antiferromagnet, while they assume op- 
posite signs on the two degenerate ground states in the 
spiral phase. In a phase with all the symmetries of the 
Hamiltonian. (Si -83) ~ (S2 -84) so that (cTa) vanishes; a 
phase with Ising ncmatic order is signaled by a {a a) ^ 0. 

Indeed, as clearly shown in Fig. 0] (lower panel), the 
critical behavior signaled by the divergence of the spe- 
cific heat corresponds to a continuous phase transition 
between a low T phase with a finite value of tr in the 
thermodynamic limit, and a homogeneous high T phase. 
Such a transition is also clearly evidenced by a sharp 
divergence of the susceptibility of the Ising nematic or- 
der parameter defined by x<t = [M /T) {{a"^) ~ 
(see Fig. 01 upper panel), and also by the T dependence 
of Binder's fourth cumulant [/4 = 1 — (ct'*)/3(ct^)^ (not 
shown) . 

The critical exponent v can be easily estimated from 
the size dependence of the T corresponding to the max- 
imum of the susceptibility, which is expected to scale as 
Tmax{L) = Tc + aL^^/'^ , where Tc is the thermodynamic 
critical temperature. As shown in Fig. 01 (upper inset), 
no sizable deviation from the Ising exponent = 1 is 
observed. The exponent /3 is also in agreement with the 
Ising universality class. This can be extracted by per- 
forming a size scaling analysis of the order parameter 
aroimd the critical T. In fact, according to the scaling 
hypothesis, close to Tc, \a\ — L~^/'^f{x), where f{x) is 
the scaling function, and x ~ tL^^^ with r = |T — Tc|. 
We have therefore plotted aL^/'^ as a function of x, using 
the value of the critical temperature Tc ~ 0.303(1) that 
can be estimated from the position of the maximum of 
the susceptibility (Fig.^ upper inset), and the behavior 
of the Binder's cumulant. As shown in Fig. 0| (lower in- 
set), an excellent collapse of the data for different lattice 
sizes on the same curve is obtained for P/v = 1/8. 

We have repeated a similar analysis for several values 
of J3 / Ji and the complete phase diagram is shown in 
Fig|Sl where we have plotted Tc vs J3/J1. We find that 
Tc vanishes linearly for J3/J1 1/4; a theory for this 
behavior will now be presented. 

Near the classical Lifshitz point, we can model quan- 
tum and thermal fiuctuations by a continuum unit vec- 
tor field n(r, r), where r = {x,y) is spatial co-ordinate, 
T is imaginary time, and = 1 at all r, t. This 
field is proportional the Neel order parameter with Sj cx 
(— l)^J+^m(rj, t). Spiral order will therefore appear 
as sinusoidal dependence of n on r. The action for n 
is the conventional 0(3) non- linear sigma model, ex- 
panded to include quartic gradient terms {h = ks = 




FIG. 5: Critical temperature as a function of the frustration 
ratio J3/J1. 



lattice spacing = 1): 5„ = j^^'^ dr J d'^rC^ with 
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— idrn) +- 
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(3) 



(dxii ■ d^n) + {dyii ■ dyii) + A2 {d^n ■ dyii 



where the ellipses denote a finite number of additional 
couplings involving 4 powers of n and 4 spatial derivatives 
invariant under spin rotations and lattice symmetries. In 
the limit ^ 00 we have x± — l/(8>/i), p ~ {Ji — 



Ah)S^, Ci = (16 J3 - ^1)^712, C2 = 0, and all A, 



0. 



Notice that p crosses zero at the Lifshitz point and so 
can be regarded as the tuning parameter; p = generally 
locates the Lifshitz point for when p < it is energetically 
advantageous to have a r-dependent spiral in n. 

A convenient analysis of the properties of 5n is pro- 
vided by a direct generalization of the 1 /N expansion of 
Rcf. [13 . The results quoted in Fig ^ and its caption 
were obtained from the N — 00 saddle point equation, 
and (apart from certain pre-exponential factors) all func- 
tional forms are exact. The saddle point implements the 
constraint = 1 and takes the form 



47r2 



Xn{k,UJ„) = 1 



(4) 



where fc is a wavevector, w„ is a Matsubara frequency, 
and Xn is the dynamic staggered spin susceptibility with 



(5) 



X„(fc,w„) = \^m^ + x±i^l + pikl + k^) 
-fCi(fc' + fc^) + 2C2fc^fc; 



The parameter m is determined by solving Eq. 1^)1. 

In the classical limit , S 00, we need only retain the 
ujn = term in Eq. 10} [1J| . A solution for m exists for all 
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T > 0, and leads to the crossovers in the spin correlation 
length ^spin shown in Fig^ The value of ^spin, and the 
pitch of the spiral order ~ \j—p^ as T — *■ are obtained 
from the spatial Fourier transform of Xn(fcjO). 

To investigate the Ising nematic order, we need to 
study correlations of the order parameter (T(r, t) which 
we define by a gradient expansion of Eq. jSJ 

a — w- dxdyti — dxn ■ dyti. (6) 

The Ising susceptibility, Xo- is then Xu = 

In the classical limit, — > cx), important exact prop- 
erties of Xt follow from the ultraviolet finiteness of 
the two-dimensional field theory with Boltzmann weight 
exp (— (1/T) / (PrCn) and n independent of r. Under a 
length rescaling analysis of this theory in which the Q 
and Xi are fixed, we see that both T and p scale as in- 
verse length squared. These scaling dimensions establish 
that in the classical limit 

6pi„ = v/cTr-SiWr) ; X. = cri$2(p/r) , (7) 

where are cut-off independent scaling function which 
depends only on ratios of the Q and A,; . The Ising phase 
transition is associated with a divergence of $2 at some 
negative argument of order unity, and Eq. Q then im- 
plies the Tc ^ —p dependence shown in Fig^ and verified 
numerically in FiglSl We can also compute Xa (including 
the quantum LOn 7^ modes) in the large N limit and 
obtain 

X. = 24r^ / ^klklxi{k,LUn) . (8) 

Using the results in Eq. (jSJ, Eq. (jSj) predicts an expo- 
nential divergence in 1/T as T for p < 0. This is, 
of course, an artifact of the large TV limit, as our Monte 
Carlo studies clearly show that Xa diverges with a power- 
law Ising exponent at a Tc > 0. 

We turn now to a discussion of the quantum physics at 
finite S. A key feature again emerges from an analysis of 
Eqs. (0J) and jSJ, while retaining the full frequency sum- 
mation: the soft spin spectrum (w ^ fc^) at the Lifshitz 
point implies that there cannot be long-range magnetic 
order over a finite regime of parameters for all finite 5* . 
After evaluating the frequency integral at T = 0, a solu- 
tion with m real exists for a range of values of | J3 — i Ji| 
smaller than ~ e^*^"^, implying there is a spin gap in this 
regime. We can reasonably expect that the Ising nematic 
order survives into at least a portion this spin gap phase, 
as it does at T > 0. 

A more careful analysis of the spin gap phase requires 
consideration of Berry phases P, Q, which are absent 
in Cn. Assuming second order quantum critical points, 
with increasing J3, we first expect a spin gap state with 
valence bond solid (VBS) order and confined S* = 1/2 



spinon excitations after leaving the coUinear Neel state. 
Conversely, decreasing J3 from the spiral spin ordered 
phase, we expect a Z2 spin liquid with Ising nematic 
order and deconfined bosonic spinous 0, 0] . So quite re- 
markably, we expect the following sequence of 4 phases 
to appear for all half-odd-integer spin S with increasing 
J3: Neel LR0-VBS~Z2 spin hquid-spiral LRO. The 2 
intermediate phases have a spin gap, and they appear 
in a window which is exponentially small in S for large 
S; the latter 2 phases have Ising nematic order. The- 
ories for the 3 quantum critical points between these 4 
phases appear in Refs. Of course, we cannot rule 

out the possibility that the some of these critical points 
and intermediate phases are pre-empted by a first order 
transition. 

It is interesting to note that other Z2 spin liquids with 
fermionic 5 = 1/2 spinous have been proposed in 
which the ground state does not have Ising nematic order. 
Our present results naturally suggest a spin gap state 
with Ising nematic order, and mean field theories for such 
states have only been obtained with bosonic spinous 01. 
Further studies of Ising nematic order in quantum spin 
models will therefore be valuable in resolving the nature 
of the spin liquid state. 
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